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The real-space renormalization group technique is introduced to evaluate the effective 
diffusion constant for diffusion in inhomogeneous media, which has been obtained by 
singular perturbation methods. Our method is formulated on a discretized real space and 
hence it can be easily combined with numerical studies for partial differential equations. 
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§1. Introduction 

The renormalization group (RG) is one of general techniques to study macroscopic (or self- 
similar) and universal characteristics of equilibrium and non-equilibrium systems whose microscopic 
models have been known. ^ 

In the present paper, the RG technique is applied to the problem of diffusion in inhomogeneous 
media. It can be considered as one of pedagogical examples of the RG method. Furthermore, since 
our real-space RG method is formulated on a discrete space, our solution can provide a simple 
example to show that the RG method and numerical simulations of partial differential equations 
can be easily combined. 

For example, there are some attempts to apply the RG technique to Navier-Stokes equation. 2 ) 
It is well known that the scale of turbulent flows of any practical significance are expected to be 
much larger than the Kolmogorov scale, which is the smallest scale of activated eddies. Hence, in 
order to save computer resource, only some large-scale eddies are computed explicitly and effects 
from other smaller eddies are modeled as a eddy viscosity. 3 ) The numerical simulations based on 
the idea of sub- grid modeling are called the large eddy simulations. The eddy viscosity used in 
large eddy simulations for turbulent fluids is expected to be evaluated with the RG applied to 
Navier-Stokes equation. 4 ) 

The only difference between the applications of the RG to Navier-Stokes equation and the 
present application to diffusion is whether the equations to be solved for decimation are non-linear 
or linear. So, our solution can be expected to help investigate applicability of the RG technique to 
solve other partial differential equations with aid of numerical computation. 

* mituhiro@riam.kyushu-u.ac.jp 
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This introductory section is concluded with note that the problem has been solved with other 
singular perturbation methods like multi-scale methods. 5 ) 

The present paper is organized as follows: After the problem to be studied is explained in § 2, 
the RG is formulated in § 3 step by step. At first, the diffusion equation is extended formally to 
close the renormalization transform in § 3.1. The construction of the renormalization transform is 
divided into two procedures, decimation explained in § 3.2 and rescaling in § 3.3. The formulae of 
renormalization transform are obtained in § 3.4. With the formulae, the fixed points are evaluated 
in § 3.5. Finally, the effective diffusion equation is obtained in § 4. The results are summarized 
and discussed in § 5. 

§2. Formulation of the Problem 

A one-dimensional inhomogeneous medium where the diffusion constant D(x) is a spatially 
periodic function is assumed. The period is denoted as I and hence D(x + l) = D(x). The objective 
of our study is evaluation of the effective diffusion constant D e when the spatial resolution of 
observation of diffusion process is much larger than the period I. The scale of observation is 
denoted as Lq. The diffusion equation is written as 

ik p{x > t] = iL D{x) iL p{x ' t); D{x +i)= D{x) - (i) 

In order to make formulation of the real-space RG easier, we discretize the continuum space 
into a chain. 8 ) The final results for the effective diffusion constant are obtained in the continuum 
limit in § 4. 

Grid points on the chain are numbered with the integers {• • • , i — 1, i, i + 1, ■ ■ ■ } and the lattice 
constant is assumed to be a (-C I), i.e., x = lim a ^oia. There are Ni sites in one period I, i.e., 
Ni = 1 1 a. The probability Pj(t) that a random walker is found at site i, or the number of the 
random walkers, is related to the probability density p(x, t) as 

p(x,t) = ]imJ5®. (2) 



Then, the diffusion equation (eq. (1)) is discretized as 

= ^ [A-iP-i(t) - (A-i + A)P(t) + DiP i+1 (t)} . (3) 

§3. Renormalization Group 

In this section, we formulate the renormalization transform for the discretized diffusion equation 
(eq. (3)) and obtain the fixed point. The construction of the renormalization transform is explained 
by being divided into two steps, decimation and rescaling. 
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3.1 Formal extension of the diffusion equation 

Before construction of the renormalization transform, new parameters are introduced in the 
discretized diffusion equation (the master equation) in order to obtain the closed renormalization 
transform as follows: 



dPi(t) 



di = h I dT [Di - l{t ~ r)p *- i(r) - - r) + A(t - r) 

+ Ui(t - T))Pi(r) + Di(t - r)P i+ i(r)] . (4) 



Although the memory effect is introduced, the effect does not exist in fact. Hence, the memory 
function Di(t) is defined with the delta function as 

Di(t) = 5(t)Di where f dr8{t - t)P;(t) = P^t). (5) 

Jo 

Furthermore, since Ui(t) is introduced only for closing the renormalization transform, Ui{t) = 0. 

In order to make decimation procedure easier, the integro-differential equation, eq. (4), is 
Laplace transformed as 

8Pi(a)-Pi(0) = ^[A-i(s)P-i(s)-(A-i(s) + A(s) + ^(s))P(s) 

+D i (8)P i+1 ( 8 )]. (6) 

Here, we denote the Laplace transform of Pi(t), Di(t), Ui(t) as Pi(s), -Dj(s), Ui(s) respectively. The 
Laplace variable s is explicitly written for prevention of confusing. For simpler notation, eq. (6) is 
rewritten as 

sPi(s)-Pi(0) = [^_i(s)P_i(s) - (wi_i(a) + Wi(s) + Vi(s))Pi(s) 

+w i (s)P i+1 (s)}, (7) 

where the Laplace transformed jump rates are defined as 

Wi{s) = Di{s)/a 2 , Vi {s) = Ui(s)/a 2 . (8) 

3.2 Decimation 

Decimation procedure, i.e., reduction of the degrees of freedom, is performed in this subsection. 
It is the first step of the renormalization transform. 

Especially in spin systems on lattices, the procedure is often implemented as averaging the 
variables in small blocks which the whole system is divided into. 9 ) On the other hand, for diffusion 
phenomena, development for long time is regarded as a kind of average process for distribution of 
the random walkers. Hence, we expect that the explicit average is not needed and reduce the degrees 
of freedom with algebraic elimination explained below. It means that the effective development 
equation for the non-averaged probability density p(x, t), instead of the averaged density, is obtained 
at last. 
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Although it appears that probability is annihilated with decimation, the renormalization trans- 
form for the initial condition, eq. (32), shows that the total probability (or the total amount of 
matter) is conserved in the renormalization transform. 

Decimation process is performed by just eliminating the probability of the nearest neighbor 
sites, Pi-i(s) and P; + i(s), from the development equation for the probability Pi(s), eq. (7). Deci- 
mation process as well as rescaling is explained schematically in Fig. 1. 

original lattice 

- a 



i-2 i-1 i i + 1 i + 2 

decimation 



i-2 i i+2 

■ *• 

i 

*» » ** 

rescaling »♦ ■ 

■ 

i jj I C i_ 

i'-2 i'-l i' i'+l i'+2 



Fig. 1. Schematic illustration of decimation and rescaling procedure performed on the discrete chain. 



Elimination of Pj_i(s) is done by solving eq. (7) where the site index is changed as i — ► i — 1. 
Pi-i(s) is given as 

p._ l(a ) = Wj-^Pi-^ + Wi-^Pijs) | Pi-i(O) ^ 

S + Wi- 2 (s) + Wi-l(s) + Vi-l(s) S + Wi- 2 (s) + Wi-l(s) + Vi-l(s)' 

The probability Pj+i(s) is obtained in the same way. By inserting Pj_i(s) and Pj+i(s) in eq. (7), 
the decimated master equation is obtained as 

Wi-l(s) 



sPi(s) - Pi(0) 



S + Wi- 2 (s) + Ifj-l(s) + Vi-l(s 



-P-i(O) 



s + u>i(s) + Wj+l(s) + v i+ i(s) 

Wi- 2 (s)Wi-i(s) 

S + Wi- 2 (s) + Wj-l(s) + Vi-l(s] 

Wi- 2 (s)Wj-l(s) 

S + Wi- 2 (s) +Wi-l(s) +Vi-l(s) 

Wi(s)w i+ i(s) 



4i(0) 

-Pi-2(S) 



S + U)i(s) + I0i+l(s) + v i+ i(s) 

swi-i(s) + 



+Vi(s) + 

4 



s + Wi- 2 (s) + ifj-i(s) + Ui-i(s) 

SWi(s) + U7i(s)Ui + l(s) 



s + u>i(s) + io i+ i(s) + fi+l(s) 

Wj(s)w i+ i(s) 

s + Wi(s) + io i+ i(s) 4- fi+i(s) 



Pi(s) 



(10) 
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Only the probabilities of the nearest neighbor sites appear in the master equation (eq. (7)). 
This is the reason why decimation is performed easily and it is locality of the renormalization 
transform. 

3.3 Rescaling 

Next, we perform rescaling procedure, which is the second step of the renormalization trans- 
form. Since one of every two neighboring sites is eliminated, the sparse grid points are renamed 
as 

»—>*' = t/2, (11) 

where i is assumed even. In the continuum limit, the rescaling eq. (11) corresponds to x — > x' = x/2, 
since the lattice constant a is not rescaled. 

Since the spatial scale is rescaled, time development slows down accordingly. It is expressed as 

i -» t! = 2~% (12) 

where the exponent \i is a positive number to be determined later. It means that the Laplace 
variable is rescaled as 

s^s' = 2^s. (13) 

Properly speaking, the probabilities {Pj(i)} are also to be rescaled. However, the diffusion 
equation that we consider is linear and hence it does not matter whether the probabilities are 
rescaled. For simplicity, we do not rescale the quantities; 

Pi(t) - pv(o = m). (14) 

Then, by using rescaling of the Laplace variable (eq. (13)) and the definition of the Laplace trans- 
form, the Laplace transformed probability Pi(s) is rescaled as 

P i ( s )^P' i ,(s') = 2-»P i (8). (15) 

3.4 The renormalization transform 

The renormalization transformed jump rates and the initial condition, w'ii(s'),v'ii(s'),P'i'(0), 
are derived by using the rescaling formulae obtained above. By inserting the rescaling formulae 
(eqs. (11, 13, 15)) in the decimated master equation (eq. (10)), the renormalized master equation 
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is given as 



s 'pV(s')--PV(o) 



2^ 



Wi^ 2 (s)Wi-l(s) 



S + Wi- 2 (s) + Wi-l(s) + fi-l(s) 
Wi- 2 (s)Wi-l(s) 



P' V -i{s') 



- T 

+2^ 



S + Wi- 2 (s) + I0i_l(s) + Vi-l(s) 

Wi(s)w i+ i(s) 



S + Wi(s) + + Ui+l(s) 

SWj-l(s) + w ._ 1 (g) 1 ,._ 1 (g) 

S + Wi- 2 (s) + lOi-l(s) + Uj-l(s) 



+2^ ( Ui(s) + 
4 



Wi(s)w i+1 (s) 



+2^ 



(16) 



s + + + v i+ i(s)' 

Since the jump rates are the coefficients of the probability P'i'(s') in the master equation, the 
renormalized jump rates w'i>(s'),v'i>(s') are given as 

Wi(s)w i+ i(s) 



w'As') = 



2 V 



S + Wi(s) + w i+ i(s) + v i+1 (s) ' 

SWi-l(s) + Wi-l(s)Vi-l(s) 



2nv i (s) + 



+ 



S + Wi- 2 (s) +Wi-l(s) +Vi-l(s) 

swi(s) + Wi(s)v i+1 (s) 



S + Wi(s) + ifj+i(s) + v i+1 (s) 

These formulae are the renormalization transform for the jump rates. 
In the renormalized master equation, eq. (16), P'j/(0) is defined as 

Wi-i(s) 



(17) 



(18) 



pV(o) = Pi(o) + 
+ 



s + u;i_2(s) + Wi-i(s) + Vi-i(s] 

Wi(s) 



Pi-i(O) 



:^+i(0). 



(19) 



s + Wj(s) +u; i+ i(s) +Uj + i(s)' 

This is the renormalization transform for the initial condition. It is weird that the renormalized 
initial condition, which should be real, depends on the Laplace variable s. However, it is shown in 
eq. (31) that the s-dependence vanishes in the limit I/Lq — > 0. Furthermore, for complete rigorous 
deduction, a new parameter should be introduced to absorb the s-dependence. If the parameter 
was introduced, the s-dependence of the renormalized initial condition would be eliminated and 
the newly introduced parameter would become zero in the limit. 

By remembering that we are interested in the behavior whose spatial scale Lq is much larger 
than /, the renormalization transform formulae obtained above are simplified. The characteristic 
time scale which corresponds to the spatial scale Lq is given by the transit time that it takes for 
the random walker to pass the region of length Lq. The transit time is of the order of Lq 2 /Di. 
Converted into the scale of the Laplace regime, s ~ -Dj/Lq 2 . From the consideration, the order of 
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s/wi(s) is estimated as 



a 

S 



a 2 D i: 



a 



(20) 



Wi(s) Df DiL 2 
Since L 3> a, s/wi(s) is a small quantity, s/wi(s) <C 1. 

The estimation of the order is not valid when the renormalization transform is performed 
many times. Here, we show that s/wi(s) is small even at the fixed point, i.e., s*/w*i* (s*) <C 1. 
Afterwards, the parameter jjl is determined as fx = 1 and hence s' = 2s. After the renormalization 
transform is performed log 2 Ni times, the fixed point is reached. Hence, the order is estimated 
as s*/w*i*(s*) ~ 0(al/Ll). Since I < L , s*/w*i*(s*) is small. It means that sW/wW i(n) (s^), 
which is the quantity renormalized n times, is small for arbitrary n, since the rescaled Laplace 
variable grows monotonically. 

The result of the order estimation means that the jump rates Wi(s),Vi(s) can be expand in 
terms of relatively small parameter s. The jump rates are expanded up to the first order O(s); 



Wi(s) = Wi + ris + 0(s 2 ), Vi(s) = Vis + 0(s 2 ). 



(21) 



The ratio of the neglected terms to the terms considered explicitly is of order of (a/Lo) 2 and hence 
contribution from these terms can vanish at the continuum limit. Furthermore, the neglected terms 
at the fixed point are relatively of order of (I/Lq) 2 . It means that the second order quantities can 
be eliminated even at the fixed point, since Lq 3> /. In addition, note that the initial values for 
Wi, ri,Vi are given by Wi = Di/a 2 ,ri = Vi = 0. 

By inserting the expansion (eq. (21)) in the renormalization transform for iOj(s), v «(s) (eqs. (17, 
18)), the renormalization transform formulae for Wi,Vi,Vi are derived. At first, the transform for 
Wi(s) is expanded in terms of s as 



w'i'(s') 

2 M 



Wi(s)w i+1 (s) 



2 M 



s + Wi(s) + w i+ i(s) +v i+ i(s) 

WiW i+ i r i+ iw 2 + riW i+ i 2 - WiW i+ i(l + v i+ i) 



Wi + W i+1 (Wi + W i+ l) 2 

On the other hand, w'i'(s') is expanded with eq. (21) in terms of s' as 

w'i'(s') = w'i> + s'r'ii = w'i> + 2^sr'i>. 
By comparison, the renormalization transform formulae for Wi,ri are given as 

WiWi+l 



+ 0(s 2 ). 



w'i> 



= 2^- 



r it = 



Wi + w i+ i 
r i+ iWj 2 + rjWj+i 2 - WjW i+1 (l + v i+ i) 

(Wi + W i+ l) 2 



(22) 
(23) 

(24) 
(25) 
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Next, by inserting the expansion, eq. (21), in the transform for Vi(s), eq. (18), the renormal- 
ization transform for V{ is given as 



«V(s')2-" 

H 4- «i, ,\ in- H 4-11, , , M 

s + 0(s 2 ) 



I0i_l(l 10,(1 
H ; h 



= u / i ,a'2-' 4 

= v'i>s. (26) 
By comparison, the renormalization transform of v; L is shown to be 

V i> = Vi H ; 1 ; . (2f) 

The renormalization transform for the initial condition is also simplified in the similar way. 
With eq. (21), the coefficients in the transform (eq. (19)) is expanded as 

Wi(s) Wi rjWj+i - Wj(l + r i+ i + v i+1 ) 



S + Wi(s) +W i+ l(s) +V i+ l(s) Wi + W i+ l (w i +w i+ i)' 2 

+0(s 2 ). (28) 

We consider the case that the transform has been performed n times. Since s ~ D^/Lq, v^ n \( n ) ~ 2™, 
the second term on the right hand side of eq. (28) is estimated as 

rW i( n)wW i( n) +1 - >Wl +r'") iW+1 +^) i( „ )+1 ) Jn) ^ yW iW j n) 



(WW i( n) +W( n ) i(n)+1 ) 2 ' wW i{ n) 



a 2 u;( n ) i( n) L 2 , 
2"aV") i(n) 

2 2n a 2 
~ -72- 



By neglecting the quantities of 0((a/Lo) 2 ),0(al/LQ),0((l/Lo) 2 ), the approximation 

^(^(S) ^(n) 



(30) 



s (n) + . (n) ( S ) + t«W i( „) + 1 (s) + »W, W+1 (s) ^( n ) i( „) + l«( n ) i( „ )+1 

is justified. Hence, the renormalization transform for the initial condition eq. (19) is simplified as 

P(" +1 W0) ~ ?W iW (0)+ (n) ""f^ P (n) ,(.)_i(0) 

+ „ t( ?\ P {n) iW +1 (0). (31) 
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Even though some variables for probability Pi(s) in decimation procedure are eliminated, it is 
shown from eq. (31) that the total probability is conserved as 

£ P(« +1 ) 1( „ +1) (0) = ^P(") IW (0). (32) 

j(n+l) j(n) 

Here, we consider the renormalized initial condition. From eq. (31), the renormalization for 
the initial distribution localized at the site miVj where m is arbitrary integer is given as 

^ (n) j(«)(°) = S i(n) jmNl / 2 n- (33) 

The result is equivalent to the limiting form at the limit t — > of the rescaling of Pi{t) (eq. (14)). 

However, the result of the transform for general initial distribution is inconsistent with the 
limit t — > of rescaling of the probability. In order to solve the inconsistency, we would had to 
reformulate decimation procedure with averaging of the neighboring probabilities. In the rest of 
the present paper, we restrict our interest in the localized initial distribution like eq. (33) for its 
simplicity. In the case, there is not any inconsistency. 

Since all the formulae are obtained, the value of \x is determined. We consider condition for 
existence of the non-trivial fixed point w* for the jump rate Wi. From eq. (24), assumption that 
the fixed point w* exists is written as 

w * = 2»~ 1 w*. (34) 

For w* is non-zero, it is necessary and sufficient that \x = 1. 

Since the number of parameters is increased, we summarize definition and the renormalization 
transform formulae: 

• The master equation 

sPi(s)-Pi(0) = [w i - 1 (s)P^ 1 (s)-(w l - 1 (s) + w l (s) + v t (s))P t (s) 

+w t {s)P l + 1 {s)} , (35) 

where Wi(s) = wi + r^s, vi{s) = vis, Di = a 2 Wi. 

• The renormalization transform 

s' = 2s (t' = t/2) (36) 

i' = i/2 {x' = x/2) (37) 

P> t ,( S ') = P l (s)/2 (P> i ,(t') = P i (t)) (38) 

w i> = 2 (39) 

Wi + W i+ l 

, _ r i+ iWi 2 + rjWj+i 2 - WjW i+ i(l + ^+i) . . 

{wi + w i+1 y 

v' v =v i+ + + Wi{Vi+1 + 1} (41) 

Wi-2 + Wi-l Wi + W i+ l 

P'AO) = P(0) + ^ Pi-i(O) + P i+ i(0). (42) 

Wi-2 + Wi-l Wi + Wi+l 
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3. 5 Evaluation of the fixed point 

At first, we consider the fixed point of the jump rate Wi. Equation (39) is converted into 
another form as 

w'i> 2 \Wi w i+ ij 

It means that the renormalization transform is equivalent to averaging of the inverse of the neigh- 
boring jump rates. Hence, the global stability of the fixed point is obvious. One period of the 
diffusion constant with Ni sites is renormalized into one grid point and the fixed point is reached, 
when the renormalization transform is iterated n\ = log 2 Z/a = log 2 N\ times. At the fixed point, 
the medium is uniform and the fixed point value of the jump rate, w*, is independent of the site. 
The value is given by 

Next, we calculate the value of V{ transformed n\ times, v*. The list of the renormalization 

transform of Vi for one period is given as 

, wi(vi + 1) w 2 (v 3 + 1) 
v i = v 2 H h 



WN, +Wl W2+ W3 



V i> = Vi-\ ■ h 



v'i' + l = V i+2 + 



Wi-2 + Wi-l Wi + Wi+l 

W i+ i{v i+ i + 1) W i+2 (Vi+3 + 1) 
Wi + W i+ l W i+2 + W i+3 

j , WjV ; -l(^jV ; -l + 1) w Nl ( Vl + l) 

v Ni/2 ~ v Nl -\ ■ 1 ■ . (45) 

W Nl - 2 +W Nl -i W N[ +W! 

Summing up the both sides respectively, the conservation law that the sum of V{ + 1 over one period 
does not change by the renormalization transform is obtained; 

Ni/ 2 Ni 

£(i/V + l) = X> + l). ( 46 ) 
v=\ i=i 

Application of the conservation law ni times gives 

2 Ni ^ 

v* + i= («&-i! + i) = -- = E(«i + i) = Ei = ^i. w 

j(«(-i) = i i=i i=i 

where the fact that the initial value for is given as v i = is used. The fixed point value is given 
as 

v* = N t - 1. (48) 

There is another parameter r; L . After iteration of the transform with n\ times, the value of 
the parameter, r*, is independent of location of the grid point. Hence, the jump rate Wi(s) is 
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transformed into w*(s*) = w* +r*s*. We show below that the second term can be neglected. Since 
s* = Nis,s/w* ~ (a/L ) 2 from eqs. (20, 36), the order of r*s* is estimated as 



\ 2 ; 
a \ „ { a 



r * Nl „ r * (49) 



Furthermore, the transform formula of and the fixed point of Vi (eqs. (40, 48)) show that 

r* ~ - A (50) 



Hence, 

(51) 



w* \Lo / 

It means that r*s* can be neglected compared to w* when the macroscopic limit I/Lq — > is taken. 

We discuss the initial condition renormalized n; times. Here, for simplicity, we assume that 
the initial distribution is localized at the origin of the chain; -Pj(O) = (5^0 • As a special case of eq. 
(33), the fixed point is given as 

P*A0) = P {ni) M) (0) = 6 i ,, o . (52) 
From eqs. (36, 37, 38), the values of i,s,Pi(s) rescaled ni times are given as 

i* = i/N t 
s* = N lS 

P*i*{s*) = Pi(8)/Ni (53) 

§4. The Effective Diffusion Equation 

The master equation at the fixed point is written as 

(v* + l)s*PV(s*) - PV(0) ^ w* [PV_i(a*) - 2P*r( 8 *) + P*i* +1 (8*)] , (54) 

where the symbol ~ denotes equality when the limits a/Lo — > 0,1/Lq — > are taken. The former 
limit is the continuum limit and the latter corresponds to the large spatial scale and the long time 
limit D(x)t/l 2 -> oo. 

With the fixed point values of eqs. (48, 52, 53), the master equation is rewritten in terms of 
the non-renormalized parameters as 

w 

N lS Pi(s) - Pi(0) ~ — [Pi-iv,^) - 2Pi(s) + P l+N[ ] . (55) 
The inverse Laplace transform gives 



dPi(t) w* 

— [Pi_ Nl (t) - 2P t (t) + P l+Nl (t)} (56) 



dt j^[P i - Nl (t)-2P i (t)+P i+Nl (t)\ 



2 * 

aw 
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The continuum limit is taken by dividing the both sides by a and the limit a — ► is taken. By 
using the definition p(x,s) = lim^o Pi(s)/a (eq. (2)) in the procedure, 

dp(x, t) D e 



[p(x -l,t)- 2p(x, t) + p(x + /, t)} 



8t ~ I 2 

is obtained. Here, we introduced the new parameter 

D e = lim a 2 w* . 



(57) 



(58) 



Since spatial scale of observation is of order of Lq, the spatial coordinate x is changed to x = x/Lq 
as 

dp(x, t) _ D e p(x - I/Lq, t) - 2p(x, t) + p(x + l/L ,t) 



dt L 2 (I/Lq) 
The limit I/Lq — > is taken with the observation scale Lq is fixed at a finite value. The result of 
the limit is given as 

(60) 



dp(x,t) _ D e d 2 p(x,t) 



dt L 2 dx 2 ' 

where we assume that the distribution of the random walker after long-time development is 
smoothed away and the derivative on the right hand side exists. Since Lq is finite and the spatial 
variable can be changed back to x = Lqx, the effective diffusion equation is obtained as 



9p(x,t) = d 2 p(x,t) 
dt e dx 2 



(61) 



It is important to note that the solution of the effective diffusion equation, p(x,t), is the non- 
averaged density. It means that the solution to the original diffusion equation (eq. (1)) at time 
t » l 2 /D(x) satisfies the effective diffusion equation without any data-processing like averaging. 

Finally, we evaluate the effective diffusion constant defined by eq. (44). The continuum limit 



gives 



D e = 



lim a 2 w* 

lim — 1 - - 



= 1/ 



[/ 

Jo 



dx 



D(x) 



(62) 



The formula, eq. (62), implies that the effective diffusion constant is mainly determined by 
regions where the value of the diffusion constant is small. In other words, the speed of the diffusion 
process is limited by the region where diffusion is slow. It can be called the bottle-neck effect. 



§5. Discussion and Summary 

In this paper, the effective diffusion constant of diffusion process in periodic inhomogeneous 
media is evaluated analytically with the real-space RG method. 
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As pointed out in § 1, the problem has been solved with the multi-scale method. We compare 
the two methods. 

In the multi-scale method, different space-time scales x, t; ex, e 2 t; . . . characterized by arbitrary 
small number e are introduced and the perturbation expansion in terms of e is performed. The 
parameter e describes separation of the two characteristic scales, the scale of the microscopic struc- 
ture of the medium and the observation scale. Hence the parameter e corresponds to the parameter 
Lq used to characterize the scale of observation in the RG method. The development equations of 
the slow variations are obtained as conditions for absence of singular and secular terms, which are 
called the solvability conditions and compatibility conditions. 

In general, the proper slow variables, e t; . . ., cannot be selected automatically and 

selection needs trial and error. On the other hand, in the RG methods the slow variables are 
obtained automatically as rescaled variables x*,t* at the fixed points. 10 ) However, the extension 
as explained in § 3.1 and the decimation procedure have to be performed in a proper way so that 
the renormalization transform is as simple as possible. Selection of the best way of extension and 
decimation needs some trial and error. 

Furthermore, the renormalization transform is formulated as recursion formulae of parameters. 
Hence, even if the transform is so complicated that evaluating the fixed point analytically is difficult, 
it would be easier to compute the fixed points numerically. Especially, the real-space RG on a 
discrete lattice space can be converted to numerical algorithms easily. 

Although the RG method has some advantages, for the problem studied in the present paper, 
calculation with the multi-scale method is much easier than that with the RG method. 
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